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Abstract 

We address the problem of the bosonization of finite fermionic systems with 
two different approaches. 

First we work in the path integral formalism, showing how a truly bosonic 
effective action can be derived from a generic fermionic one with a quar- 
tic interaction. We then apply our scheme to the pairing hamiltonian in 
the degenerate case proving that, in this instance, several of the features 
characterizing the spontaneous breaking of the global gauge symmetry U(l) 
occurring in the infinite system persist in the finite system as well. 
Accordingly we interpret the excitations associated with the addition and 
removal of pairs of fermions as a quasi- Goldstone boson and the excitations 
corresponding to the breaking of a pair (seniority one states in the language 
of the pairing hamiltonian) as Higgs modes. 

Second, we face the more involved problem of a non-degenerate single particle 
spectrum, where one more kind of excitations arises, corresponding to the 
promotion of pairs to higher levels. This we do by solving directly the 
Richardson equations. From this analysis the existence emerges of critical 
values of the coupling constant, which signal the transition between two 
regimes, one dominated by the mean field physics, the other by the pairing 
interaction. 

1 Introduction 

The problem of the bosonization of finite and infinite fermionic systems has been 
extensively explored in the past. All these investigations have been prompted by 
the recognition that in nature this phenomenon is widely occurring. Best examples 
of it are offered by the superconductivity in metals with the related BCS theory 
and by the superfluidity of certain atomic nuclei [2] where it is signalled by the 
existence of a gap (about 1 MeV) in the energy spectrum, by the reduction of about 
a factor of two in the moment of inertia with respect to the shell model prediction 
and by the staggering behavior of the separation energies (odd-even effect in the 
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mass number A). More generally in the nuclear case the bosonization phenomenon 
is of course epitomized by the outstandingly successful Interacting Boson Model 
(IBM) by Arima and lachello P|. In this paper we first address the problem of the 
bosonization in the path integral formalism. Although this approach is general and 
hence keeps its validity in bosonizing any fermionic action with a quartic interaction, 
for sake of illustration we shall consider the well-known pairing hamiltonian which 
reads 

L ju L 

v=l mu=—ji/ fi,u=l 

In O 

CLj^m^ is the destruction operator of a fermion in a single particle level (whose 
number is L) characterized by angular momentum j^, third component and 
single particle energy 6^, whereas the operator 

Ai, = ("I)"'" "^^ O'j^-m^^j^m^ (2) 

m„=l/2 

destroys a pair of fermions with total angular momentum J = in the level jiy. 
The Hamiltonian is indeed: 

1. relevant for both nuclear and condensed matter physics. In nuclear physics it 
is known to represent, together with the quadrupole force, a dominant part of 
the residual nucleon-nucleon interaction and a crucial microscopic ingredient 
of the IBM i3j. If the interacting fermions are electrons, H (in particular Hp) 
is the Hamiltonian underlying the BCS theory of superconductivity ^ H|. 
Moreover, the same Hamiltonian governs collective phenomena in the physics 
of liquids and metal clusters; 

2. a simple example to illustrate the bosonization of fermionic systems. The 
problem of relating the fermionic Hamiltonian to a bosonic one, like the 
IBM, with the same spectrum, has been the object of several investigations [5| 
ini ISj and still represents a challenging question. 

Furthermore: 

3. for infinite systems, the BCS theory, of which the hamiltonian offers a 
simple realization, is an example of spontaneous symmetry breaking of the 
gauge group U (1), with the associated appearance of a Goldstone boson jjj, set 
up with excitations consisting of the addition and removal of pairs of particles 
to and from the system. The question then arises whether this Goldstone mode 
survives in finite systems and, if so, how the Goldstone and Higgs excitations 
should be identified. 

4. The Hamiltonian has been shown to be exactly integrable [S], which means 
that the number of constants of motion equals the number of degrees of free- 
dom. This does not imply, however, that explicit expressions for its eigenvalues 



Bosonization of the Pairing Hamiltonian 



3 



and eigenstates are easily obtained, since the constants of motion are usually 
complicated operators which, in general, can only be diagonalized numerically. 
For this reason it is useful to provide approximate analytic expressions for the 
eigenvalues and eigenvectors of Hp. 

This report is organized as follows. In Section |2l we illustrate how the path 
integrals scheme works and apply it to in the degenerate case, where only one 
single particle level is considered, deducing the well-known zero seniority spectrum 
in the path integral formalism, discussing the appearance of a Goldstone boson and 
studying the seniority excitations within a matrix approach in Sec. 12.51 In Section |31 
we study, by directly solving the Richardson equations, the problem of one and 
two pairs living in many non-degenerate single particle levels, deriving analytic 
expressions for the energies of both the collective and the trapped states. Finally, 
in Sec. |l]we briefly discuss the relation between the exact and BCS solution to the 
pairing problem in terms of Bogolioubov quasi-particles. 



2 The degenerate case 

The pairing Hamiltonian has been first diagonalized by Kerman, Lawson and Mac- 
Farlane [H] in the simple case of one degenerate single particle level by using group 
theoretical methods. In fact H can be cast in the following form 

L L 

H = 2j2euSzu + Yj^'^^'^-GS+S- , (3) 

u=l u=l 

= jv + 1/2 being the pair degeneracy of the level ji,, in terms of the quasi-spin 
operators piU ITT] 

L L 

5+ = J] it, S. = Y,^u, (4) 

v=l u=l 
^ L L 

u=l u=lmv=l/2 

which span a SU(2) algebra. Hence the pure pairing Hamiltonian Hp can be imme- 
diately diagonalized, yielding 

-G < S+S- >= -G [N/2 {n - N/2 + 1) - v/2 (0 - v/2 + 1)] , (6) 

where N is the eigenvalue of the particle number operator (the number of fermions 
is assumed here to be an even integer) 
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and = the total degeneracy. The seniority quantum number 

v = Q-2S , (8) 

S{S + 1) being the eigenvalue of S*^, represents the number of unpaired particles. 
Now, in the degenerate case {L = 1) = e and the full Hamiltonian becomes 

H = 2e {S, + n)- GS+S^ . (9) 

Hence the spectrum generating algebra is still SU(2) and the energy turns out to 
read 

E{n, s) = 2en-G[n{n-n + l) - s{n- s + l)] . (10) 

In (jlOj) . for future convenience, the pair number n = N/2 < Q and the quantum 
number 

s = v/2 < min{n, Q — n} , (11) 

which counts the number of broken pairs, have been introduced. 

The spectrum (fTU)) is associated with two independent types of excitations: one 
is related to the addition or removal of one pair of fermions and the other to the 
breaking of a pair. The former, described by the quantum number n, is a Gold- 
stone boson associated with the spontaneous breaking of the global gauge invariance 
reflecting the particle number conservation; the latter, described by the quantum 
number s, can be viewed as corresponding to the Higgs excitations |12| . 

In an infinite system the energy of the Goldstone bosons vanishes with the asso- 
ciated quantum number. This does not occur in a finite system, but, to the extent 
that the energy spectrum of the latter displays a pattern similar to that of an infinite 
system, it should exhibit two quite different energy scales. Actually the excitation 
energies associated with both the quantum numbers n and s appear to be of the 
same order, namely gQ. However the Goldstone nature of the energy spectrum as- 
sociated with n is clearly apparent when one considers the excitations with respect 
to the minimum of pOj) . This occurs for 

n = no = M, (12) 

[...] meaning integer part, with 

^o = ^(^^ + l)-^. (13) 

Introducing then the shifted quantum number 

u = n — no , (14) 

(fTU]) becomes 

E{no + u,s) = Gu"^ + 2Gu{no - uq) - Gno{2uo - uq) + Gs{n - s + 1) . (15) 

Therefore the addition (or removal) of one pair of nucleons with respect to the 
ground state requires an energy of order G: this is the energy of the Goldstone 
boson. The energy required to break a pair, namely the seniority energy, is instead 
of order Gil: this is the energy of a Higgs boson. 
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2.1 The path integral approach 

In this Section we show that the pairing spectrum for s = can be obtained in 
the framework of the Feynman path integral and illustrate how an effective bosonic 
action can be set up in this case. 

To do this it is necessary to introduce the odd Grassmann variables Am, AJ^ 
associated to the fermionic destruction and creation operator am, aln (we omit the 
index j since we are considering only one single particle level) and the even elements 

^mit) = (-1)^ — A_m(t)Am(t), ^*mit) = i-iy~"'X*mm-mit) , (16) 

describing pairs of fermions with = 0. 

Since we restrict ourselves to the zero seniority part of the spectrum ()1U() . we 
define 

m = i-iy-^X.mmmit) , (17) 

m>0 

namely the even elements of the Grassmann algebra associated to the J = com- 
posites (jSJ, also referred to as "hard bosons". The analysis of seniority excitations, 
which will be studied in 12.51 in a different framework, would clearly require the 
introduction of J 7^ composites. 

Note that the index of nilpotency of the composites ()17|) is 0, (i.e., <I>" = 0, 
Vn > fl), reflecting the Pauli principle. 

In terms of the above variables the euclidean action associated to the Hamilto- 
nian reads 

No/2-1 f \ 
S = ^ E \Y.>^Ut){yt + e)Xrn{t-l)-G^''m{t-l)\, (18) 

t=-NQ/2 I m ) 

where the time has been discretized (r is the time spacing), A'^o is the number of 
sites of the lattice, the discrete time derivative 

vfm = ±-[f{t±i)-f(.t)] (19) 

r 

and e the single particle fermion energy. The variables Am at the time t = —Nq/2 
are related to the ones at t = No/2 — 1 by antiperiodic boundary conditions 

Am(-A^o/2) = -Am(iVo/2 - 1) . (20) 
The action H18|) can be written in terms of the even variables as 

No/2-1 f \ 

= E E [-VlXt)Vm{t) - X^ip*^{t)^m{t - 1)] - GT^*m{t - 1) , (21) 

t=-No/2 lm>0 J 
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where 

X = 1 - re . (22) 

Following the approach of Ref. for performing Berezin integrals over com- 
posite variables, the generating functional and the correlation functions of n spin 
zero pairs of nucleons can then be expressed as integrals over the (f variables of the 
composites, according to 

Z = j [dip*d^]e'^ (23) 

and 

< <I>"(t2) [<^*{h)r >=\j Wd^W{t2) WihTe-^ , (24) 
respectively. In the continuum limit the following result is then obtained 

< $"(t2) [^*(ti)]" > ~ e-^/^^-^i)^" (25) 

for the correlation function, with /3i.2 = Tti^2 and 

En = n[2e - G{n - n + 1)] (26) 

for the energy spectrum. The latter exactly coincides with ()1U|) when s = 0: hence 
the well-known zero-seniority formula for the spectrum of the system is recovered. 

Concerning the bosonization of the fermionic system, although in the pairing 
Hamiltonian the fermionic and bosonic degrees of freedom get mixed, a purely 
bosonic effective action for the field $ can be deduced using the variables above 
introduced. It reads IT] 



t=-No/2 lfc=l J 

(27) 

ak being numerical coefficients whose asymptotic behavior is found to be 



Oik 



. (28) 



The action 5*^*1 is equivalent to the fermionic action (|18j) in the sense that it leads 
to the same spectrum ()1U() . Worth noticing is that its structure is characterized by 
the absence of terms of the type [^* {t)f[^{t)f-P[^t - l)Y (p / 0), due to the 
occurrence of nontrivial cancellations. 



2.2 Goldstone bosons in the Hubbard-Stratonovitch approach 

In what follows we again consider the part of the spectrum associated with the 
addition and removal of pairs, showing its connection with a Goldstone boson. 
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It is well-known that the spontaneous breaking of the local electromagnetic gauge 
invariance plays a crucial role in the theory of superconductivity jl5l 171 ITB]. Here 
the symmetry group U(l) of the global gauge transformations ip e^^ip is broken 
by the non- vanishing expectation value of pair operators (the Cooper pairs), which 
carry charge —2: the symmetry group is thus reduced to the unbroken subgroup Z2, 
consisting of the two gauge transformations with a = and vr. The Goldstone field 
associated to this symmetry breaking lives in the coset space U(l)/Z2 and displays 
only derivative interactions. 

In a finite system, like the one considered here, a spontaneous symmetry breaking 
cannot rigorously occur, since tunnelling takes place between the various possibly 
degenerate states and the true ground state turns out to be a unique linear super- 
position of the degenerate states 0. However, if the two above mentioned features 
of the Goldstone fields in infinite systems, namely that they parameterize the coset 
space U(l)/Z2 and have only derivative interactions, survive in finite systems, then 
the identification of a Goldstone boson becomes of relevance not only for a deeper 
understanding of the bosonization mechanism, but also for a convenient choice of 
the variables. 

To explore this occurrence we perform a Hubbard-Stratonovitch linearization of 
the action H18() . introducing an integration over auxiliary fields in the generating 
functional 

Z = J [d\d\*dr]dr]*]e-^ , (29) 

where the new action is 

No/2-1 r 
t=-No/2 I 

+ E K{t){Vt + e)Xrn{t)] \ (30) 

m=-j J 

and the auxiliary bosonic fields rj* and r] satisfy periodic boundary conditions. 
Next we introduce the following polar representation for the auxiliary fields 

r? = ^e'^', n* = V~Pe~^'' (31) 

and explore whether the field can be identified with the Goldstone boson. 

First, for the change of variable (|3H) to be one to one (with the only exception 
of the point p = 0), 9 must vary in the range < < vr: hence the field 6 lives in 
the coset space of the broken symmetry group U{\) of particle conservation with 
respect to the unbroken subgroup Z2, as appropriate to a Goldstone field [7]. 

Second, the Goldstone field should display only derivative couplings in the action, 
which is not the case for the field 9 after the transformation (|31j) . However the non- 
derivative coupling can be eliminated introducing the following transformation on 
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the nucleon fields: 

\ — je j \* _ I* 
In fact, as a consequence of the above, the fohowing operators 

= exp(=Fi^)V^ exp(±i6') ± e , 
whose time matrix elements are 



exp jir (V^+e)^^ } St^^t^+i - 6tit2 
5tit2 - exp jir {Vf 0) 5t^ ^t^^i 



- e5tit2 , 



(32) 
(33) 

(34) 
(35) 



will appear in the action. In the latter the 9 field appears only in these operators 
and therefore under derivative, as appropriate to a Goldstone field. 

2.3 The saddle point expansion 

Having proved the Goldstone nature of the excitations associated to the addition or 
removal of fermionic pairs, we now proceed to deduce the energy spectrum through 
the saddle point expansion. 

After integrating over the fermionic fields ipm and ip^ we get for the generating 
functional the expression 



(36) 



with 



[dp] / [d9]exp{-S,ff), 



Seff = r^Gp-Trln{-q-q+ + G^p) , 



(37) 



where the trace is meant to be taken over the quantum number m > and the time. 
The U{1) symmetry is now realized in the invariance of Seff under the substitution 



+ (38) 

with a time independent. 

To perform the saddle point expansion one should first look for a minimum of 
the effective action ()37() at constant fields. Denoting with p the time-independent 
component of the p field and defining 

M = + G^p (39) 



and 



p-i = -(l-er) V+V,-+M2, 



(40) 
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we can write 5'e// at constant fields as follows 

5e// = r^Gp-TrlnP-i . (41) 
t 

By performing the trace over m and converting the sum into a time integral, taking 
first the Nq oo limit and then letting r ^ (with tNq constant), the above 
becomes ^21 

5e// = tNq {Gp + ne- nU) , (42) 
whose minimum with respect to p occurs when 

M = ^, (43) 

or, equivalently, at 

P = P^^ (^ [(^^)' - ^^'] = ^ ■ 

It is remarkable that the dimensionless 'pq, when multiplied by G^, coincides with 
the well-known gap A characterizing the BCS theory in the one single particle level 
case. The action Seff at the minimum is then 

So = rNo^-^-'^+ne^ . (45) 

Next we perform an expansion around the saddle point. We start by defining 
the fluctuation of the static p-&e\d according to 

P = Po + r = Po(l + ^) (46) 
and by considering the generating functional (jHSJ, now written as 

/oo rn 
[dr] / [d9]expi-Seff). (47) 
-Po -^0 

Obviously this expansion is justified only if the quantum fluctuations in (|46() are 
small. We will show at the end of this Section that this is indeed the case. 
To proceed further we rewrite S^ff in the form 

5e// = r^G(po + + TrlnP-Trln[l + P(i?i+i?2)] , (48) 
t 

where 

Ri = -q'q+ + {V+ + e){Vt -e) (49) 
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and 



We set then 



i?2 



s. 



eff 



oo 

r=0 



(50) 



(51) 



the term 5o being the saddle point contribution, given by H45() . This grows Uke 0^, 
but it contains also a term of order J7 and a term of order one, which should be kept 
if an expansion in powers of 1/Vt is sought for. However, we prefer to stick to the 
definition 1)40(1 for the operator P and to compute the further contributions to the 
expansion (|51|1 (the quantum fluctuations) by developing the last logarithm in the 
right-hand-side of (|48l) : the terms thus obtained are naturally organized in powers 
of . It is worth noticing that this expansion does not break the U{1) invariance. 

Actually we shall confine ourselves to consider the first and second order contri- 
butions in Vj0 and r. 

The first order action stems from the term linear in r and from the first term in 
the expansion of the logarithm and reads 



Si = TGY,rt-TT[P{Ri + R2)] 



(52) 



The explicit computation shows that all the terms linear in r cancel out: hence the 
r-integration remains undefined. However the contributions arising from the second 
term in the expansion of the logarithm make the integral over r well defined. The 
second order contribution will not be reported here and the details can be found in 
Ref . • One is thus left with the result 



Si = -Tv{PRi) = i 



l + (M + -e)T 



(53) 



Note that these contributions are of order fl and 1, namely are 0{1/Q) with respect 
to the saddle point one. 

In order to obtain the Goldstone boson energies we must find out how they 
depend upon the single particle energy e. For this purpose we have to perform in 
the integral expressing the generating functional Zi (associated with the action Si ) 
the ^-integration, which is not gaussian, because 6 is compact. Yet we can choose 
VtO as a new integration variable, thus rendering the integral gaussian, getting 



Not 



■InZi 



3 M 



3e + Gn 



(54) 



Finally, in order to deal with a specific system (e.g. a nucleus), the particle 
number must be fixed. This can be accomplished by replacing e with 



e=e- n , 



(55) 
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/i being the chemical potential and using 

Id Id 
<N>= — lnZ= T^lnZ, (56) 

where N is the particle number operator. Now, replacing < N > with 2n and 
noticing that M does not depend on /x (since it is independent of e, see Eq. ^^), 
Eq. becomes 

which gives 

;x = G(n-J7/2-3/8) + e (58) 

for the chemical potential. Hence, in the presence of the chemical potential, the 
energy of the system becomes 

En,o = ;^(5'o - InZi) + 2/in 

= 2en - Gn{n - n + 3/i) + - {5n + -] , (59) 



which reproduces the excitation spectrum of the pairing hamiltonian with good 
accuracy, the relative difference between the exact ground state energy —G[{Q + 
l)/2 — e/G]^ (namely Eq. ()1U() with s = and n = no) and the one corresponding 
to Eq. (|59|) being of order 

Note that the same result can be obtained, rather than through the chemical po- 
tential, by introducing in the path integral the particle number projection operator 

Vn = p ^e-^(^-2")" , (60) 

as shown in Ref. jl2j . 

An important comment is in order on the validity of our expansion, which de- 
pends upon the size of po- We assess the latter by replacing in (|^H) e by e and using 
(|58jl (dropping the irrelevant term -3/8), thus obtaining 

-pQ = n{n-n) , (61) 

which attains its maximum value pg = ^^^/4 for n = Q./2. This corresponds to the 
situation when the level where the pairs live is half-filled. When this situation is 
almost realized, namely when the shell is neither fully filled nor almost empty, the 
functional integral defining Z becomes gaussian and an expansion in r/pQ can be 
performed. On the other hand po attains its lowest value when n = 1 or n = Q. In 
this case there is no shift in (|46j) and the p field acts only through its fluctuations, 
which are small, thus supporting the validity of the expansion. It is then not sur- 
prising that our approach yields the (almost) correct excitation spectrum for any 
value of the pair number. 
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2.4 The hamiltonian of the s-bosons 

We are now in a position of deriving the bosonic hamiltonian corresponding to the 
effective action previously obtained. The most general, particle conserving, quartic 
hamiltonian for a system of s-bosons, confined to live in one single particle level, 
reads in normal form (we omit for brevity the indices specifying the single particle 
levels on which the b operators act) 

H{b\b) = hiy^b + h2bWbb, (62) 

b\ b being bosonic creation-annihilation operators acting in a Fock space and satis- 
fying canonical commutation relations. 

The parameters /ii,/i2 should be fixed by requiring that the spectrum of (|62() 
coincides with the pairing Hamiltonian one. This is indeed possible and it turns out 
that the choice 

hi = 2e-Gn , h2 = G (63) 

accomplishes the job. However the Hamiltonian thus obtained, being intrinsically 
bosonic, patently violates the Pauli principle and therefore the condition n < 
should be added a posteriori, when using (|62|) in dealing with a system of fermions. 

On the other hand in our framework this condition naturally emerges. Indeed if 
we write the path integral associated to (|62|) 

Z = J [db*db]exp{-S), (64) 

where 

No 

s = Tj^[b:+iVtbt+Hib:^„bt)] , (65) 

t=l 

b*,b being holomorphic variables satisfying periodic boundary conditions in time, 
and perform the same saddle point expansion previously illustrated, from the com- 
parison with our effective action we get 

hi = 2e-GQ- G/4, h2 = G, n< n . (66) 

The inequality in the above equation, expressing the Pauli principle, is necessary for 
the two path integrals (|36() and (|64() to coincide and follows from the positivity of po 
as given m Eq. ((nU: thus in our approach this condition, far from being artificial, is 
necessarily implied by the formalism itself. Obviously the considerations following 
(|5H) hold valid here as well, hence hi will be affected by an error of order 1/M. 
Accordingly the values (|66|1 essentially coincide with (|63j) . thus giving the exact 
spectrum. 

Summarizing, we have constructed an effective bosonic Hamiltonian, namely 



H{b\b) = (2e - Gn)b^b + Gb^b'^bb, 



(67) 
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which yields the exact zero-seniority pairing spectrum and naturally incorporates, 
when treated within the path integral approach, the Pauli condition n < ^l, keeping 
track of the fermionic nature of the composites. 

2.5 Seniority excitations: a matrix approach 

In this Section we address the problem of the whole spectrum of the degenerate 
pairing Hamiltonian, including the excited states with non zero seniority v. For 
convenience we set e = and consider only the interaction Hamiltonian Hp. 

The path integral approach followed in the v = case has not been pursued till 
now for V ^ since it becomes quite cumbersome. Hence in this case it is preferable 
to work in the hamiltonian formalism. 

In the framework of the creation and annihilation operators the commutator 

i, if 

is non-canonical, thus rendering not trivial to find the eigenstates of Hp. An ap- 
proach circumventing this difficulty consists in employing again even Grassmann 
variables: in so doing the difficulty associated with the non-canonical nature of the 
commutator (|68() disappears, since the associated Grassmann variables do commute. 

Within this formalism we study whether and how the Fock basis, set up to 
with determinants of single particle states, can be reduced to a minimal dimension 
without loosing the physical information we search for and whether this minimal 
basis can be expressed in terms of composite bosons [T7| ITS]. 

In seeking for the reduction of the basis dimensions, we shall be guided by the 
two main features of Hp, namely that it is expressed solely in terms of the (p and is 
invariant for any permutation of the ip. 

These properties, indeed, urge us to express also the vector of the basis in terms 
of the variables (|16j) . The action of Hp on these states is then obtained following 
the lines illustrated in ^1], the result being 

HpiPiip*) = J [dip'*dip']Kp{ip\ip')e^'^'*'^'iP{ip'*) = EiPiif*) , (69) 

where the integration is over the even elements of the Grassmann algebra and the 
kernel reads 

Kp{y,*,ip') = Hp{p\^')e^^*^' . (70) 

We then attempt to diagonalize the Hp associated with n pairs in a basis set 
up with states represented as products of n factors ip*^s. Since the number of these 
is (^) , the very large reduction of the basis dimension entailed by the choice of the 
variables (|16j) is fully apparent: indeed in terms of fermionic degrees of freedom, the 
corresponding basis would have had a dimension (g^) . 
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Moreover, since the variables A anticommute, each vector of this basis is anti- 
symmetric in the exchange of any pair of fermions, thus fulfiUing the Pauh principle. 

We now actually explore whether, for a given and n, the eigenstates of Hp 
can be cast into the form of a superposition of products of n variables ip*, namely 

o 

V'(<^*) = E/^™['^-i---'^™J™' (71) 

m=l 

the index m identifying the set of quantum numbers {mi,m2 • • • m-n} find the Prn 
being complex coefficients. 

From the condition 1)11(1 it follows that the number of distinct eigenvalues of Hp 
is given by [min{n, — n} + 1], each of them having a degeneracy ^H] 

s being, we recall, the number of broken pairs (in our convention a binomial coeffi- 
cient with a negative lower index vanishes). 

We look for the eigenvalues of Hp by diagonalizing the symmetric matrix Tip = 
(E — Hp)/G, hence dimensionless, given by 

(£+n OVl 
OVl S+n 

The dimension of the above is (^) and the symbol OVl indicates that the upper 
(lower) triangle of the matrix is filled with zeros and ones. Indeed the matrix 
elements oiUp are one when the bra and the ket differ by the quantum state of one 
(out of n) pair, otherwise they vanish. The diagonal matrix elements simply count 
the number of pairs and £ = E/G. 

We now look for a further dimensional reduction of the basis such that the 
resulting matrix Hamiltonian has the same eigenvalues of the original one, however 
all being non degenerate. For this purpose an elementary combinatorial analysis 
shows that the number of ones in each row (column) of the matrix ()73() is given by 
n(0 — n). Indeed a non- vanishing matrix element has the row specified by n indices 
whereas, of the indices identifying the column, n — 1 should be extracted from those 
fixing the row in all the possible ways, which amounts to n possibilities. The missing 
index should then be selected among the remaining Vl — n ones: hence the formula 
n(il — n) follows. 

Note that the cases with n and 17 — n pairs are equivalent, an occurrence which 
is also manifest in the exact spectrum (|1()|) . Hence in the following we shall confine 
ourselves to consider n <^ only. 
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Figure 1: The figure shows, in the specific case n=3 and 17=9, the partition of the 
O levels in the set I (with 17 — n=6 levels) and II (with n = 3 levels). 



To write down explicitly the matrix (|7l-ij) it is convenient to divide the set of 
the even Grassmann variables, whose quantum numbers identify the levels where 
the n pairs are placed, into two subsets: one with — n and the other with n 
elements (to be referred to as I and II, respectively). This partition leads to a basis 
with dimension n + 1, independently of the value of fi. The associated eigenvalues 
correspond to the breaking of 0, 1, 2 • • • n pairs. 

Indeed in this instance the (^) entries of each row and column of the matrix can 
then be grouped into n + 1 sets, the first one corresponding to the n pairs placed in 
the — n levels of I, the remaining n levels of II being empty (see Fig. la). 

The number of configurations belonging to this first set is do = {^~^)- To the 
second set are associated configurations with n — 1 pairs in I and one pair in II (see 
Fig. lb), their number being di = ''^(^Z")- In general the {k + l)-th set embodies 
configurations with n — k pairs in I and k pairs in II, their number being 

4 = ^ 2^ with < A; < n (A; integer). (74) 
Clearly the total number of configurations is 



5^4 

fc=0 



n 



(75) 



With this organization of the levels the matrix (|73j) splits into (n + 1)^ rectan- 
gular blocks Bkj (with < k,j < n) of dimension dk x dj and, since the pairing 
Hamiltonian only connects states differing by the quantum number of one pair, the 
blocks with \k — j\ > 2 will have vanishing elements: the matrix thus becomes 
block-tridiagonal. 
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This suggests the introduction of the following orthonormal set of n + 1 even 
commuting variables (the composites) 

^i = yr £ [<-'^-Jm ' (76) 

where k varies as in ((7H) . Sj = X^^^q being s_i = and m again identifies the set 
{mi, 1712 ■ ■ ■ TTT'n}- Notc that the variables (|76|) have in general an index of nilpotency 
higher than one. 

The definition (|76() also reflects our desire that the composite bosons keep as 
much as possible of the symmetry of Hp. And indeed the while not fully 
symmetric with respect to the interchange of the ip*, turn out to be invariant with 
respect to the interchange of the (/?* belonging either to the set I or to the set II. 
This symmetry property enforces the maximum coherence among the components 
of It is remarkable that composite variables corresponding to combinations of 
the ip^_^ (^^2 ■ ■ ■ Vtrin different from ((75|) not only hold a lower symmetry than the 
one displayed by (fTHj) . but may also lead, as we have verified in some instances, to 
the wrong eigenvalues. 

Now in the minimal basis H76|) Tip is represented by a (re + 1) x (n + 1) matrix 
whose generic element (M„)^^ =< ^'^\Ti.p\^* > obtains by summing the elements 
of the block B^i of the (^) x (^) matrix, but for the normalization factor l/y/dyh- 
The sum is performed by recognizing that all the blocks have the same number of 
ones in each row. Specifically, in the diagonal block B^k the number of "ones" in 
each row is (re — k){^ — 2n + 2k). In the upper diagonal block Bk^k+i each row 
contains instead 

Cfc = (re - kf (77) 

ones. Since the total number of ones in each row of the matrix H73() is re(r2 — n), the 
number of ones in the rows of the lower diagonal block B^^k-i will be 

bf, = k{n -2n + k) . (78) 

As a consequence the non vanishing elements of the matrix M„ turn out to be 



{Mn\k+i = <'^k\np\n+i>=^j^^c,, (79) 

(M„),+i_, = <^l^^\np\n>=^^b,^i (80) 

and 

(.Mn)kk =< ni'^pin >=ak = £ + n + n{n-n)-bk-Ck . (81) 



Bosonization of the Pairing Hamiltonian 



17 



Clearly (M„)^^_,_-^ = (M„)^_,_j^ ^, since the operator Tip is Hermitian and the basis 
H76|) is orthonormal. The matrix thus becomes tridiagonal, reading 



M„ 



\ 



ao 



ai 





02 





ak 



C2 



dfe+1 







-1 an J 



and the associated eigenfunctions should be expanded in terms of the composite 
variables (|76jl . namely 



(82) 



fc=0 



It is convenient to cast the coefficients of the expansion (|82() , fixed by the eigen- 
value equation 



MnU = , 



into the form 



dkWk 



(83) 



(84) 



Thus the (^) equations for the /5's reduce to n + 1 equations for the w^s and the 
associated eigenvalues £ obey the secular equation 



Dn = det (M„) = , 



(85) 



where £ enters into the diagonal matrix elements afc. 

Now from the general theory of symmetric tridiagonal matrices one knows that 
Eq. H85() has n + 1 distinct and real roots and these are found by applying the 
recursive relation 



Dn = an.Dn-l (Cn-1 )^-Dn-2 



(86) 



for increasing values of n. Hence we have 
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Do = y 

Di = y{y-n) 

D2 = y{y-n)[y-2{n-l)] 

Dn = y{y-^)[y-2{^-l)][y-^{a-2)]...[y-n{^-n + l)] , (87) 

where 

y = n{n-n)+£ + n . (88) 

We thus see that Dn has all the zeros of -Dn-i plus an extra one for y = n(J7 — n + 1). 

Moreover 1)8 7|) allows us to write down for the general solution of 1)85(1 the ex- 
pression 

y = p(n -p+1) with <p<n . (89) 
Hence the well-known formula for the spectrum of the pairing Hamiltonian 



£ = -{n-p){n-n-p+l) , (90) 

is recovered, the index p coinciding with the pair seniority quantum number s. 

Let us now consider the eigenfunctions of Hp. For a tridiagonal matrix, a 
recursive relation among the eigenvectors components, similar to (|86() . can also be 
established. In the specific case of the matrix M„ it reads 

dkCkWk+i = -dk-iCk-iWk-i - dkakWk , (91) 

where again < k < n and quantities with negative indices are meant to be zero. 
Thus for the lowest eigenvalue y = {p = 0, zero seniority) we have 

ujt"^=wt'^ = .... = u:t'K (92) 
namely the collective state 

o 



^.^0 = p E Vd-,n = Q E [^*m.-^Um ■ (93) 

k — TTh — 1_ 



m. 



Indeed in ()93|) all the components of the wave-function, i.e. the monomials [(p'^^...(p: 
are coherently summed up. This state obtains for a specific partition of the levels 
defining the matrix M„. However, any other partition would lead to the same result, 
being all the weights of the components equal. As a consequence the state (jHSJ is 
non degenerate. 



nJ m' 
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Concerning the second eigenvalue y = $7 (p = 1, pair seniority s = 1), according 
to Eq. (|91() the components of its eigenstate are 



w 



,^ -J\fi{kn-n^) . (94) 

Finally the components of the state associated with a generic pair seniority s turn 
out to read 

j=0 \j) 

= Ms 3F2{-k,-n + s-l,-s] -n,-n; 1) (95) 

3F2 being a generalized hypergeometric function. In the above the binomial (J) is 
meant to vanish when j > k. 

In particular for the vector of the maximum seniority w^^^"^^ , corresponding to 
y = n(0 — n + 1) {s = n), one has 



{n-2n + k)\{n-k)\ 
^ ' {n-2n)\ 



(96) 



In and (|^^ A/i, Ms and Mn are normalization constants. 

It is interesting to reobtain in the present formalism the following important 
result j2Uj : for values of large with respect to n, the number of dominant compo- 
nents in w^^'^ decreases with s, reflecting the weakening of collectivity with increasing 
seniority. Indeed in the limit 17 >> n the eigenvalues are 

y ~ sl^ (97) 

and moreover 

Cfc << 57 , Cfc„i ~ kQ and ~ (s — fc) 17 , (98) 

being < s < n. 

By inserting the above limits in H91|) the w^^^ components corresponding to the 
eigenvalues (|97j) are found to be 

t(;Q = = ... = w)._i = (99) 

and ^ 

4'^ = -y-^l'^ fc = 5,...,n. (100) 



We thus see that the state with seniority s has indeed, in the basis of the ydfe^fc 
and in the large $7 limit, s vanishing components, the remaining n — s + 1 ones being 
expressed, through (|100|) . via the single component Wn^. In other words, when 17 is 
large, the collectivity of a state with seniority s decreases as s increases, because its 
components become fewer and fewer and, furthermore, the surviving components 
are more and more expressible through a single one. 
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3 The non degenerate case 

We now turn to explore the case of the pairing Hamiltonian H acting on a set of 
non degenerate single particle levels. The actual situation one faces in applying the 
pairing Hamiltonian to real systems corresponds indeed to having more than one 
energy in In this respect nuclei and metals represent two extreme situations 
of the non-degenerate case: in the former a major shell is typically split into five 
or six single particle levels of different angular momenta, in the latter the number 
of non-degenerate levels entering into a band corresponds to a significant fraction 
of the Avogadro number. Moreover in a heavy nucleus the number of pairs living 
in a level may be as large as, say, eight while in a metal is one. Recently, a re- 
newed and widespread interest for the pairing problem in the non-degenerate frame 
has flourished in connection with the physics of ultrasmall metallic grains, possibly 
superconducting 21 , and of Bose-Einstein condensation |22j . 

As we shall illustrate, in the non-degenerate case one more type of excitations 
occurs with respect to the degenerate one, corresponding to promoting pairs above 
the Fermi sea: these are zero seniority excitations which, in the BCS language, 
are described as four quasi-particle states (more precisely, as two quasi-particles 
and two quasi-holes states). These are conveniently classified in terms of a new 
quantum number, the "like-seniority" , and reflect the occurrence of a quantum phase 
transition. 

Furthermore, and more importantly, in the non-degenerate case critical values 
Gcr of the coupling arise which split the physics into a region governed by the mean 
field (when G < Gcr) and a region governed by the pairing force (when G > Gcr)- 
This competition does not show up in the degenerate case. 

3.1 The Richardson exact equations 

In this Section we derive in an alternative way the well-known Richardson equa- 
tions |23j . using the Hamiltonian formalism expressed in the framework of even 
Grassmann variables previously introduced. As already mentioned, this is partic- 
ularly useful for avoiding the non canonical commutators, which would naturally 
appear in the standard Hamiltonian treatment of the pairing problem. 

In the following we consider the simpler case of zero-seniority states, although 
the approach can be generalized to include seniority excitations by introducing other 
composite variables, as done for one pair in Ref. j24j . 

Starting from the normal kernel of the Hamiltonian written in terms of 
Grassmann variables 




(101) 



we search for eigenstates of n pairs of fermions in the s-quasibosons subspace as 
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products of n factors, namely 

n 

M^lim) = llBl{m) , (102) 

k=l 

where 

L 

is a superposition of s-quasibosons placed in all the available levels and the index 
m = {mi, . . . , TJin) labels the unperturbed configuration from where the state devel- 
ops as the pairing force is switched on. The set of values of m corresponds to the 
possible states available to the system. When no confusion arises the index (m) will 
be dropped. 

It is convenient to start from the effective Hamiltonian 

L ji. L 



coincident with (|101|) in the s-quasibosons subspace spanned by the states (|102j) . 
Indeed while terms like A*A count the number of particles, ip*ip = A*A*AA counts 
the number of pairs. The eigenvalue equation can then be written as 

j [V*dV^']Wefr((^*,</^')exp f ^^(V.*^^^ +</.;*^J(p;.^„J Vn(^'*) = i^nV'n(^*) , 



(105) 

since in the expansion of the exponentials only the even powers, hence only the if 
variables, survive. By performing the integrals over the (p's one gets then 

where Ck{m) are normalization factors (see for their expression) and Ek{m) 
(henceforth referred to as "pair energies") are the solutions of the non- linear system 

y._iV y i ^1 k = l,.--n. (107) 

^2en-Ek ^ Ei-Ek G ^ ' 

We thus see that the Grassmann variables formalism enables us to recover the same 
equations obtained within, for instance, the quasi-spin framework and usually known 
as Richardson's equation. 

Since the equations (|1()7|) deal with pairs of fermions (e.g. nucleons) coupled to 
an angular momentum J = 0, their eigenvalues, given hy E = Ej., are those of 
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the zero-seniority states {v = 0). Importantly, these eigenvalues display different 
degrees of coUectivity: hence they are conveniently classified in terms of the latter. 

To clarify this point we introduce a number vi, that counts in a given state 
the number of particles prevented to take part into the collectivity, not because 
they are blind to the pairing interaction (indeed they are coupled to J = 0), but 
because they remain trapped in between the single particle levels, even in the strong 
coupling regime. This number is directly linked to the number Nq, first introduced 
by Gaudin in Ref. through the relation vi = 2Ng and might be considered as 
a sort of "like-seniority" (hence the notation vi) ^ since it reduces to the standard 
seniority v for large G: as pointed out in |^, Nq has the significance of the number 
of pair energies which remain finite as G goes to infinity. 

Specifically we shall ascribe the value vi= to the fully collective state, vi=2 
to a state set up with a trapped pair energy while the others display a collective 
behavior, f;=4 to the state with two trapped pair energies and so on. 

Worth noticing is that in the degenerate case e^ = e the result ()1U|) for vanishing 
seniority is easily recovered by multiplying Eq. p()7|) by {2e — Ek) and summing over 
k. Indeed, exploiting the identity 



which coincides with (jlOj) for s = 0. 

In the non degenerate case an exact analytic solution of the Richardson system 
cannot be obtained and one has to resort to numerical calculations. However, ap- 
proximate analytic expressions for both the ground and the excited states can be 
obtained and they help in understanding the occurrence of critical phenomena. This 
will be illustrated in the next two Sections for the cases of 1 and 2 pairs. 

3.2 One pair 

We start by addressing the problem of one pair living in many single particle levels 
and show that, even in this critical value of the coupling constant G can 

occur which separates two different regimes |28j . 

For n = 1 the Richardson system (jlUTjl reduces to the equation 



^Note however that, whereas the seniority v counts the fermions coupled to J 7^ 0, the number 
vi refers to J = pairs. 




(108) 



one immediately gets 




(109) 




(110) 
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which yields L eigenvalues E{m) {1 <m < L), while the wave function 

L 



is characterized by the coefficients /3j/(m) = C{m)/{2ei, — E{m)), according to 
Eq. (ITHHl) . 

It is straightforward to solve equation (|11U|) numerically: the solutions can be 
graphically displayed as the intersections of the left-hand-side of UllUf) with the 
straight line E = 1/G PHI- Two classes of states appear: the first one embodies the 
lowest energy state, which lies below the lowest single particle level for an attractive 
interaction and represents a collective state when G is large with respect to the 
single particle energies; the other contains the so-called "trapped" solutions, which 
lie in between the single particle levels. 

Concerning the wave function, the one corresponding to the lowest eigenvalue, if 
G is large enough to organize a collective motion, has all the /3 coefficients sizable, 
reflecting the high degree of collectivity of the state. On the contrary, in the eigen- 
states of the trapped solutions only one /3-coefficient dominates, so that their wave 
functions are close to the ones of an unperturbed state. 

Let us now consider the collective state in the strong coupling regime. Here it 
is convenient to recast Eq. piOj) as 

V - I V — = - ai2) 

V u 1 — Z— 

E-2e 

and expand in the parameter 2(ei, — e)/{E — 2e). Defining 



in leading order we get 



(113) 



Eo = 2e-i}G, n = ^n^, (114) 

which coincides with the degenerate case value. Owing to the definition (|113() . the 
next-to-leading order correction vanishes. To proceed further we rewrite H112() as 

n - n Gn Y ^ i 



E - 2e ^ (E - 2e)'' E-2e'^\E-2eJ G 

n=0 n=0 ^ ' 

where the generalized moments of the distribution of single particle levels e^, namely 

Min) = ^n^in) = i ^ n,{e, - eT , (116) 
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have been introduced together with the dimensionless expansion parameter 
and the variance 

^ = iii^y:^Aeu-er . (ns) 




The strong couphng regime then corresponds to a <^ 1. Note that the second 
moment is the square of the variance a, whereas the third, M^^^ = cr^j, and fourth, 
M^^^ = (c+3)fT^, moments are related to the skewness 7 and to the kurtosis c of the 
distribution, respectively. Note also that a = aM~^, being the dimensionful 
expansion parameter introduced in Eq. (|43)1 . 

In a perturbative frame, setting E^^^ = £^("~^) + S and linearizing in 6, we get 
for the lowest energy 

E = 2e-Gn[l + a^ --ra^ + {l + c)a^ + 0{a^)] , (119) 

an expression valid for a <^ 1 and for any single particle energy distribution, but in 
particular when the fluctuation of the latter around e is small. 

The evaluation of the trapped eigenvalues can be easily performed numerically. 
This approach however hides some interesting properties of the solutions, which 
correspond to a transition from the mean field to the pairing domain, as we shall 
illustrate in the following. 

The number of trapped solutions of the secular equation (|11()|) for one pair and 
L single particle levels is L — 1. They can be identified by a quantum number v = 
2, . . . L, which labels the level where the pair lives when G = (v = 1 corresponds 
to the lowest unbounded energy). Since they are trapped, namely 

2e^_i < E^""^ < 2e^ , (120) 

a new variable z^'^^ G (0, 1) can be defined according to 

i^^'^) =2e,_i + 2z(^)(e,-e,_i) . (121) 

If we now isolate in the secular equation the terms associated with the poles in 
2eu-i and 26^, we get 

^u-i _ 2(e^ — e^,_i) noo") 

li{fv,v—l)=l z'' ' 

Next we let the discrete variable v become continuous by means of the Euler- 
McLaurin formula [30) . which replaces the sum with an integral, i.e. 



6 

E 



f{y) ^ l/(a) + \f{h) + f{u)du . (123) 
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The formula H123|) holds valid if /(i^), i' to be viewed as a complex variable, is 
analytic in the strip a < Rez^ < b {a and b being integers). 

To proceed further one must specify the single particle energies and the associ- 
ated degeneracies. Here we shall first consider the case of a harmonic oscillator well, 
and later generalize the results. The harmonic oscillator unperturbed energies are 
labeled by an index k = 0,1, ■ ■ ■ M — I and read 

ek = ekhojo = {k + 3/2)nwo , {k = 0,1, ■ ■ ■ Af - 1) (124) 

with associated pair degeneracies 

= (A; + l)(/c + 2)/2 . (125) 
The total degeneracy Q, the average energy e and the variance a are found to be 

n = lMiM+l){M + 2) , (126) 



and 



a 



e = -{Af + l)nuJo 



2- ^ {M-l){M + 3){nLOo? , 



80 



respectively, entailing the following expression for the expansion parameter: 



G V 5AA(AA + l)(AA + 2) 



(127) 
(128) 

(129) 



Using the dimensionless quantities G = G/fuvQ and z = z/hwQ we then rewrite p22|) 
in the Euler-McLaurin approximation obtaining (A: labels the single particle levels) 



n{k-i) ^ n{k) 



z{k) - 1 



E-McL 



(k,~z{k)) 



(130) 



with 



R-McL 



{k,z) 



+ 



- ^(4A: + 2z + 3) + ^(AA - 1){M + 2A: + 2z + 3) 
k(l-k) (/c + 2)(A: + 3) 



4(5 + 1) 
1 



4(z - 2) 
M{M + 1) 



2(A; + 5-l) 4(/c + z-AA) 



+ -ik + ~z){k + ~z + l)\o^ 



(z + l)(A: + i-AA) 



{z-2){k + z-l) 



(131) 
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The Euler-McLaurin approximation provides the key for studying analytically the 
limiting case — > oo, which is helpful in shedding light on the properties of the 
trapped eigenvalues. For this purpose we define 



A 



k 



(132) 



which, when J\f is large, tends to become a continuous variable in the interval (0, 1) 
Moreover we express the coupling constant G in terms of a according to H129() 
Then Eq. (|130j) becomes 

m = eoo(A) 



with 



and 



1 



1 _J_ 

^ 2-1^ 2(z + l) ^2(2 



1 1 - A 1 - 

^o.(A) = - + log^ + — 



log 



5 + 1 



z-2 



a 



(133) 
(134) 

(135) 




Figure 2: Solution of the eigenvalue equation in the limit ^ oo in the Euler- 
McLaurin approximation for different values of a. Solid line: a = 0, dashed line: 
a = 0.2, dotted line: a = 0.35, dash-dotted line: a = 0.5, long-dashed line: a = 0.8, 
dot-dot-dashed line: a = 1.5. The exact results are not displayed as they almost 
coincide with the ones in the figure. 
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The numerical solutions of Eq. ()133() are displayed in Fig. [21 for different values 
of a. It is interesting to follow the behavior with a of the curves, from the a — > 
case (strong coupling), which still carries the fingerprints of the harmonic oscillator, 
to the a — > oo one (weak coupling), which corresponds to the straight line z = 1. 

The behavior of the curves is ruled both by the function ^oo(A), which goes to 
— oo when A — > 1, whereas for A ^ 

Coo(A) '-^ I +~ « < with a„ = lS^ 0.29 , (136) 

[ — oo tor a > Ocr 8 V 5 

and by the monotonic decrease of 4>{z) in the interval (0, 1), at whose endpoints it 
assumes the values 

lim 0(5) = +00 and lim (f){z) = — oo . (137) 

Thus, when A — > 1, 

no longer depends upon a and all the curves in Fig. [2 coalesce to 1. Indeed since 
the pairing interaction is of finite range, a pair trapped in highly excited harmonic 
oscillator states has the two partners sufficiently de-localized to be little affected by 
the interaction. In this connection it is worth reminding that the classical limit is 
achieved by letting the degeneracy of the single particle levels become very large 
|31j . Thus the coincidence of all the eigenvalues, for any G, in A=l also reflects the 
evolution from quantum to classical mechanics of our system. 

For A — > two cases occur: if a < then i(A) — > 0, whereas if a > 
then z(A) — > 1. Thus a transition occurs at illustrated in Fig. |2| the almost 

parabolic behavior of z{X) for small a (which was first observed in Ref. [2^) is 
strongly distorted for a ~ Oct for small A. For larger a a smoother behavior is 
recovered. In particular in Fig. |2a marked minimum is seen to develop for a above, 
but close to, the critical value. 

The convergence of all the curves of Fig. [2 to z(0) = for a < Ocr reflects the 
pressure exercised by the infinite number of the high-lying, large degeneracy, levels 
on the low-lying, low-degeneracy, ones. This occurrence might be understood on 
the basis of a sum rule the trapped solutions should fulfill and of the nature of the 
pairing force. Actually the sum of the eigenvalues 5(A) can be exactly computed in 
the two limiting cases a = and a = oo: from the Viete equations one indeed has 

= 0) ^ E ^(^) = i (139) 

k=l 



and 



E(a = oo) = 1 , 



(140) 
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respectively. These values set the limits for the area under the curves of Fig. [2J 
The existence of a critical point reflects the fact that, since the action of the pairing 
force is gauged by the product Gfi^, above some critical value of a the system 
prefers to obey the sum rule by lifting the lowest eigenvalues (corresponding to the 
lowest degeneracies) to the unperturbed values. Indeed, as we shall illustrate, if the 
degeneracy of the single particle levels is not growing fast enough with k, then no 
transition occurs. 

The case of finite M will not be discussed here. We just mention that, as proved 
in Ref. |22], the eigenvalues obtained in the AA — > cxd are very robust with respect 
to variations of M, when G is large: indeed they keep their validity even for values 
of as small as five. 

It is now natural to ask whether the transition previously discussed is peculiar 
of the harmonic oscillator model or is more general. To answer this question, we 
take J\f large and consider two classes of models with equally spaced single particle 
energies, as in the harmonic oscillator case, but degeneracies growing (model a) or 
decreasing (model b) with k, namely 

n^^^ = (A; + 1)T ~ AA^A^ (141) 

nf^ = {AT- ky = Af^ii - \y , (142) 

with 7 > 0. The associated expansion parameters are then 



V (7 + 3) (7 + 2)GAA7 ^ 

Note that for 7=2 and M large the harmonic oscillator expression for a is recovered, 
up to a factor of two. 

By performing again the Euler-McLaurin approximation and studying the be- 
havior of the trapped solution a transition is found to occur only for the models a, 
with 

4?^ = , (144) 

" V 7 + 3 7(7 + 2) ' ^ ^ 

clearly behaving as I/7 when 7^0. On the contrary, for the models b no transition 
can take place. Thus a transition occurs only if the degeneracy grows with k, and the 
'strong coupling' domain becomes wider (ocr increases) as 7 approaches zero: here 
the transition disappears. If the degeneracy decreases with k (case b), no transition 
exists. 

3.3 Two pairs 

We now address the problem of two pairs living in many levels, which is already 
quite instructive on the general case of many pairs and on the evolution towards the 
BCS limit. 
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In the previous Section we have shown that even when the system consists of 
only one pair a transition between two different regimes is found, a precursor of the 
quantum phase transition occurring in infinite systems. When the system is made 
of several pairs, this phenomenon becomes more compHcated, since many critical 
values of the coupling G arise. These critical values correspond to a particular 
"escape" mechanism of the trapped solutions from the grid of the single particle 
levels, which is necessary for the occurrence of a fully collective state, as we shall 
now illustrate. 

As we shall see, as for the one pair case, a transiton between two different regimes 
still occurs, but now the number of critical points can be, under suitable conditions, 
two, due to the larger number of possible configurations. If the number of single 
particle levels tends to infinity the two critical values of G merge into one, which, 
notably, coincides with the G^ of a one pair system. In correspondence of this Gcr 
the system undergoes a transition from a mean field to a pairing dominated regime. 

The Richardson equations (|1U7() reduce, for n = 2, to the following system of 
two equations 

^=1 "^e^ — E2 El — E2 

As discussed in Section |31 the solutions are classified in terms of like-seniority: we 
shall ascribe the value vi= to the fully collective state, vi=2 to a state set up with 
one trapped pair energy and one displaying a collective behavior and v;=4 to the 
state with two trapped pair energies. 

Let us first consider the weak coupling limit, where of course no collective mode 
develops (hence vi has no significance). Adopting a perturbative treatment for the 
pair energies Ei, E2 we write 

Ei = 2e^, + Gxi , (146) 

Gxi being a perturbation. Then 

\^ = --^ + Y" ^ (147) 

and, expanding in G, the system (|145|) becomes 



f 

1 + 






Xi 


< 

1 + 






X2 



(148) 



+ ^ + 0{G') = 0, 



where the indices fii and ^2 select one configuration out of the unperturbed ones. 
At the lowest order in G (weak coupling regime), if ^ui 7^ fi2 (namely if the two pairs 
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sit on different single particle levels when G = 0), one has Xi = —^^^ and the pair 
energies Ei = 2e^. — Gfl^. are real. Thus the total energy E = Ei + E2 becomes 



If /ii = /i2, which is possible only if VL^-^ > 1, the Richardson system is 



1 + 







Ml 



1 + 



Xl 

Q 



X2 



n,. 



+ 



'^^'^i 2(e^ - e^J X2-XI 



n 



+ 0{G' 



and from its solution 



xi,2 = -{^fii - 1) ± iV^^n ~ 1 



(149) 



(150) 



(151) 



the complex conjugate pair energies Ei = E2 are obtained. The total energy reads 
then 

E = 4e^, - 2G{n^, - 1) (152) 

and is, of course, real. 

In comparing (|149j) and l|152j) with the energy of one pair system in the weak 
coupling regime, namely 

E = 2e^- GQf, , (153) 

one sees that, while (|149j) corresponds to the sum of two contributions like (|153j) . 
hence the two pairs ignore each other, in (|152() a positive (repulsive) energy 2G 
associated to the Pauli blocking appears. 

Let us now consider the strong coupling domain. We shall deal only with the 
states vi = and 2 since the f ; = 4 states are of minor physical interest. 

The vi = state arises from an unperturbed configuration with the two pairs in 
the lowest single particle level (if l^i > 1) or in the two lowest ones (if f^i = 1). 

Following the same procedure as in the one pair case, we introduce the variables 



Ei-2e 

Gn 



(154) 



and solve the system (|145|) in terms of the expansion parameter a defined in (|117|) 
through a recursive linearization, getting 



3^1,2 



a 



{n-r 



(^-4) 

(n-i) 



lb ia 7- 



(155) 



{n - 1)3/2 

Note that the two pair energies are always complex conjugate and the system's total 
energy reads 

2n{n - 2) 



E = 4e-2G{n - 1) 



1 + 



(n-i) 



(156) 
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Hence the leading order contribution coincides with the degenerate case collective 
eigenvalue (see Eq. [TU)) , namely 

^deg = 4e - 2G{Q - 1) , (157) 

which therefore represents a good estimate when the spreading of the mean field 
levels is small with respect to Gfi, as it is natural to expect. 

We note that, when 0^1, (|156j) becomes just twice the value (|119|) of the 
collective energy of one pair of nucleons living in L levels in the strong coupling 
regime. Moreover, the imaginary part of xi and X2 goes to zero as l/^/Tl. Thus, 
in this limit, the Pauli interaction between the two pairs vanishes, as expected: the 
two pairs behave as two free quasi-bosons condensed in a level whose energy is given 
by (tTTTHl . 

Let us now come to the like-seniority vi = 2 states. 

If the coupling term (expressing the Pauli principle) in the Richardson equations 
(|145|) were absent, then the eigenvalues of the vi = 2 states could be simply obtained 
by adding the collective energy E\ carried by one pair and the trapped energy E2 
carried by the other pair. 

This situation is recovered in the very strong coupling limit, where all the single 
particle energies become essentially equal to e and both e and E2 are very small 
with respect to Ei. The first equation of the system (|145() then becomes 

I n 2 , , 

r + Tr--Er = (158) 

yielding 

El = -G{n - 2) , (159) 

namely the energy of the state with two pairs and v = 2 in the one level problem. 
This result sets a correspondence between states with vi = 2 and v = 2, thus 
connecting seniority and "like- seniority" (or the physics of a 'broken' and a 'trapped' 
pair). 

In the strong coupling domain however the Pauli term in (|145|) cannot be ne- 
glected, but is well approximated by 2/Q: hence the system H145() decouples and 
can be recast as 
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where 



G 



1 



'(1) 



1 




2 




(161) 



cff 



1 



1 



G 



1 



- E- 



2 



(1) 




(162) 



1 



Therefore the Pauh principle in the large G regime just re-scales the coupling con- 
stant in such a way that the pairing interaction is quenched for the collective state 
and enhanced for the trapped ones. 

The results obtained in the strong coupling domain are in agreement with the 
findings of Ref. [Slj, where the problem of the pairing Hamiltonian for small super- 
conducting grains is studied. The main difference with the present approach is that 
there the single particle levels are non-degenerate: hence the results of Ref. j32j for 
two pairs can be recovered from ours by setting = I. 

Let us now discuss the existence of one or more critical values for the coupling 
constant G. 

In the weak coupling regime, when a state evolves from an unperturbed one with 
the two pairs living in the same level, then the pair energies Ei and E2 are always 
complex conjugate. On the other hand when the state evolves from an unperturbed 
one having the two pairs living in two different levels, then Ei and E2 are real. 

By contrast, in the strong coupling regime the pair energies Ei and E2 of the 
vi = state are always complex conjugate. It is then clear that, if the degeneracy 
f^i of the lowest single particle level is greater than one, then the pair energies Ei 
and E2 are complex conjugate in both the weak and strong coupling regime and 
their behavior with G is smooth. 

On the other hand if f^i = 1, since in the weak coupling limit the two pairs 
must live on different levels, Ei and E2 are necessarily real in a neighbourhood of 
the origin, but become complex in the strong coupling regime. Thus a singularity 
in their behavior as a function of the G is bound to occur. 

In this second case it appears natural to surmise that the singularity takes place 
when El and E2 coincide. In fact in this case the Pauli term of the Richardson 
equations diverges (reflecting the attempt of four particles to sit in one level that 
can host only two of them), and it must be compensated by a divergence in the sum 
entering into the system (|145|) : this can only happen if Ei or E2 coincides with an 
unperturbed eigenvalue. 

The situation is portrayed in fig. 01 where the behavior of the pair energies with 
G is displayed for a harmonic oscillator well assuming L = 3. 

Clearly in this case two values of Gqt occur. The pair energies Ei and E2, real 
in the weak coupling limit, coincide at the critical point G^r^^ (their common value 
being 2ei) and then become complex conjugate; the energy of the associated state 
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"1,2 




2 - 




0.2 0.4 0.6 0.8 1 1.2 

G 



Figure 3: Behavior with G of the pair energies of the ground {vi = 0) and first 
excited (f/ = 2) states obtained as solutions of the Richardson equations for a 
harmonic oscihator well for L = 3: dashed lines denote the common real part of 
the two pair energies when they are complex conjugate and solid lines the real part 
of the two pair energies when they are different. All energies are in units of Hujq. 
As explained in the text the two critical points refer to the = (G ~ 0.21) and 
-t;; = 2 (G ~ 0.53) states. 



evolves in the vi = collective mode. By contrast, for the vi = 2 state, arising from 
the configuration with two pairs living on the second level (which is allowed for the 
harmonic oscillator well), the two pair energies Ei and E2 are complex conjugate 
in the weak coupling limit, coalesce into the energy 2ei at the critical point 
and then become real. One of the two solutions remains trapped above 2ei , while 
the other evolves into a collective state: the sum of the two yields the energy of the 
vi = 2 state. 

Analytic expression of the critical values of G can be obtained starting from the 
Richardson equations and assuming 0.^ = 1. Indeed in Ref. j25j the expression 



i-r 



(2)1. 



E 



n, 



E 



'2e„. - 2e, 



(163) 
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where G"*" is the critical value of the couphng related to the ground state and G 
the one related to the excited states, is derived. In the above the quantities 



are the inverse moments of the level distribution. We thus recover, although through 
a somewhat different route, the results found long ago by Richardson |23j . 

It is important to stress that for the occurrence of a critical value of G a single 
particle level with $7 = 1 must exist. Such level is the lowest one in a harmonic 
oscillator: for this potential well both the ground and the first excited state of a 
n = 2 system carry one critical value of the coupling constant, namely Gcr^^ and 
gH^ , respectively. In correspondence of these values of G the pair energies take 
on the value Ei = E2 = 2ei . Thus for a n = 2 system two (at most) critical points 
exist on a = 1 level. 

Finally, owing to the relevance of the = 1 degeneracy, we consider the model 
of L, for simplicity equally spaced, single particle levels all having = 1, a situation 
occurring in metals and in deformed nuclei. In this instance two positive Gcr always 
exist in the lowest level when L > 3 (in fact G^J^ — > 00 for L = 2). Moreover 
a positive G^^ implies complex Ei and E2 for G < Gcr^ and, since for small 
G both the pair energies are real, they should evolve from an excited unperturbed 
configuration having the two pairs placed in the two lower-lying single particle levels, 
as illustrated in figlH Numerically, for this model, we have found that two Gcr appear 
on the second level when L > 9 and on the third level when L > 16. Thus in the 
Q = 1 model two pairs can form, so to speak, a quartet only if they live on adjacent 
single particle levels in the unperturbed configuration. Furthermore the more excited 
the unperturbed configuration is the more not only G, but also L, should be larger 
for the merging to occur, a fact clearly reflecting the competition between the mean 
field and the pairing force. Finally we observe that here, at variance with the finding 
of Ref. imiSIlES], Gir^^^+ > G^r^+: this is simply because we measure the single 
particle energies from the bottom of the well rather than from the Fermi surface. 

It is interesting to establish the relationship of the above discussed critical values 
of G with the critical value (|136() found for one pair, which corresponds to 



when L is large. For the n = 2 system we first notice that, at large L, in Gcr the 
moment Vi dominates over V2- hence, in this condition. 




(164) 



G, 



cr — 



_8_ 
I2 



(165) 



~ (7(1) 

^cr cr 



(166) 



Bosonization of the Pairing Hamiltonian 



35 



10 




Figure 4: Behavior of the pair energies of three excited states in the case of = 1. 



The number of levels considered is L 
of the level spacing. 



9. The energies £'1,2 are expressed in units 



Moreover, again for L large enough, 

^1 ^ Y , (167) 

entailing the equality 

Gcrin = 1) = Gcr{n = 2) (168) 

in the asymptotic L limit. From this outcome one sees that also when n = 2 the 
relevant dynamical element for G < Gcr''^ is the mean field, whereas for G > Gci^ 
it is the pairing force, as far as the system's ground state is concerned. 

3.4 Transition amplitudes 

Further insight into the critical behavior of the 2-pair system can be gained by 
studying the pair transfer matrix elements from a one pair to a two pairs state as a 
function of G. 
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This can be achieved by constructing the wave functions according to Eqs. 1)102(1 . 
(|1U3|) . (|1U6() with n = 1 and n = 2 and, with these, calculate the transition amph- 
tudes induced by the operator (see Eq. ((21)), namely the matrix elements 

<n=l\A^\n = 0> and <n = 2\A^\n=l> . (169) 

In the fully degenerate case the general expression for these transition amplitudes 
between states with any number n of pairs and seniority v, can be obtained using 
the expression Q for the excitation spectrum. Indeed, inserting into the pairing 
spectrum a complete set of states \n',v' >, one has 

^ < n + l,v\A^\n',v' >< n',v'\A\n + l,v >= + 1 " ^) (^^ - n - 1 - ^ 

n',v' 

(170) 

thus getting for the transition amplitude 

<n+l,v\A^\n,v >= + 1 " ^) " ^ " ^) • (l"^!) 

Note the G-independence of ()171|) . which, moreover, is diagonal in the seniority 
quantum number. From the previously illustrated analogy between seniority and 
like-seniority, we expect that for large G the matrix elements (|169() display an asymp- 
totic behavior coinciding with (|171|1 . namely 

< l,v, = 0|it|o > > ^/n (172) 

,1+, f a/2(J7 - 1) fort;, = , , 

G^oo \^yu — 2 toT vi = 2 . 

Actually this happens only if vi is conserved. However, the number vi is not con- 
served at finite G: only in the strong coupling regime the vanishing of the matrix 
element between states of different Gaudin numbers occurs. Indeed for finite G we 
get 



< (mi, 1712)1 A^\{m) >= 
^E[/^?^("^i'"^2)/3(^)(mi,m2)]' 



, 2 

v ^ V 



(174) 



M being a normalization constant. 

Using the Richardson equations to get rid of the sums, we recast ()174|) as follows 

M^ti^ ^^ ^ 2Ci(mi,m2)C2(mi,m2)Ci(m) 
< {mi,m2)\A^\im) >= ——— , 175) 

N G[{Ei[mi,m2) - E{m)\[{E2{mi, m2) - E{m)\ 
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Figure 5: The transition matrix elements ()175() divided by ^/Tl for L =4 (solid), 10 
(dashed), 20 (dotted) and 40 (dot-dashed) levels of a h.o. well. Upper left panel: 
transition from the m = 1 (collective) state to the (mi, 7712) = (1,2) {vi = 0) state; 
upper right panel: transition from the m = 1 (collective) state to the (mi, 7712) = 
(2, 2) {vi = 2) state; lower left panel: transition from the m = 2 (trapped) state to 
the (nil, 7712) = (1,2) (vi = 0) state; lower right panel: transition from the m = 2 
(trapped) state to the (mi,m2) = (2,2) {vi = 2) state. An asterisk denotes the 
position of the critical points. The straight lines represent the asymptotic values for 
G ^ 00. 
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which is real since Ei and E2 are either real or complex conjugate. 

We display in fig. [3] the amplitudes (\175\i for the transfer from a n = 2 to a 
n = 1 system (divided, for obvious convenience, by We consider the n = 1 

system either in the ground collective state (m = 1) or in the first excited trapped 
state (m = 2). On the other hand the n = 2 system is either in the vi = (namely 
(mi, 7712) = (1,2)) or in the vi = 2 ((mi,m2) = (2,2)) state. The calculation is 
performed for a harmonic oscillator well with L = 4, 10, 20 and 40. 

The matrix elements clearly exhibit a conspicuous enhancement in the proximity 
of Gcr (marked by an asterisk in the figure), which most clearly reflects the onset of 
an ODLRO (off-diagonal-long-range-order) in the system. 

The behavior of the transition matrix element around the critical points is quite 
delicate since both the numerator and the denominator in H174() tend to vanish when 
El and E2 tend to 2ei. In this limit the transition amplituse (|174j) reads 

<(mi,m2)|^t|(^)> ^ (176) 

G^Gt. Gcr ^6PJ) ± 8P(2)Pf3) + 2PJ) 

Now, the inverse moments of the level distribution (|164() vanish in the L — > 00 limit: 
hence (|176j) diverges. As above mentioned, this occurrence, not appearing in fig. [3 
because of the chosen normalization, relates to the ODLRO which sets in into a 
system close to a phase transition. 

Specifically the diagonal amplitudes for large L behave according to 



< (mi, m2)|A"''|(m) >'^~°° 



/217 G - Gcr = 



The off-diagonal amplitudes behave instead as 

< (mi,m2)|i"f|(m) >^~°°(5((5-Gcr) • (178) 

From the figures it is also clear that the critical value of G increases with L and 
Ga^'^ ^ 8 as L ^ 00 (see Eq. ((Tl 



4 BCS and the exact solution 

In concluding this work we shortly comment on the relation between the exact 
solution of the Richardson equations and the BCS solution in terms of Bogoliubov 
quasi-particles. 
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Figure 6: Behavior with G of the excitation energy of a two-pairs system in a 
harmonic oscillator well with L = 3. The exact solutions with vi = 2 (circles) and 
vi = 4 (triangles) are compared with the BCS results for u = 2 (solid line) and v = 4 
(dashed line). For each like-seniority the two exact solutions develop from different 
unperturbed configurations. 



To this scope we self-consistently solve the well-known BCS equations for N 
fermions living in L levels, namely 



" 2 

L 



V(e.-/x)2 + A2 



1 - u^, 



E 



L 



1 



V(e,-/x)2 + A2 



(179) 
(180) 
(181) 



where A = G^^^ O.uUuVu the gap, fi the chemical potential and = eu — Gv'^, and 
we evaluate in this framework the excitation energy of a system with seniority v, 
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given by the energy of v quasi-particles, namely 

V 

Ebcs = \/(e.-A^P + A2 . (182) 

v=l 

In fig. iniwe display and compare the exact excitation energies -Eexc = E{vi) — 
E{g.s.) for a L = 3 harmonic oscillator well and the corresponding Bogoliubov's 
quasi-particles predictions ()182() for a v = 2 and a v = 4 state. It appears that for 
G larger than the highest critical value (G^ ~ 0.6) both the BCS and the exact 
excitation energy become linear functions of G and, remarkably, are very close to 
each other, in particular for vi = 2. 

Note that the excitations in the BCS picture amount to break a pair, whereas 
in the Richardson frame the like-seniority excitations correspond, microscopically, 
to two particles-two holes (the vi = 2) and to four particles-four holes (the vi = 4) 
excitations without the breaking of any pair. In other words they are associated 
with the promotion of one or two pairs to higher lying single particle levels. 

The result in Fig. El thus strengthens the correspondence between seniority and 
"like-seniority". Moreover it confirms that the Richardson exact solutions behave 
as (|182() in the strong G limit, as was proved by Gaudin in the large L limit [^Hl- It 
is remarkable that this appears to be approximately true already for L = 3. 

5 Conclusions 

In this report we have presented two approaches to the bosonization of fermionic 
systems based on the introduction of composite variables associated with the rel- 
evant degrees of freedom. One relates to the change of variables in the Berezin 
integral expressing the generating functional, the other to the straight solution of 
the Richardson equations. 

Our testing ground has been the simple pairing Hamiltonian, which acts between 
fermions coupled to total angular momentum J = 0. Such system exhibits three 
kinds of excitations, related to the addition of removal of pairs, breaking of pairs 
(the so-called seniority excitations) and promotion of J = pairs to higher energy 
levels, respectively. 

The first kind of excitation connects the ground state energies of different sys- 
tems, in particular atomic nuclei, and is associated with a Goldstone mode in the 
spectrum. Although the system we consider is finite, and therefore cannot truly 
display a spontaneous symmetry breaking mechanism (as it happens instead in su- 
perconductivity), nevertheless the precursor signature of a Goldstone boson can be 
clearly identified. This has been demonstrated, at least in the degenerate case, by 
performing a Hubbard-Stratonovitch linearization of the action and showing that 
the associated bosonic auxiliary field indeed endows the properties of a Goldstone 
boson. 
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The seniority excitations have been studied, again in the degenerate case, in 
the framework of the Hamiltonian formahsm cast in the language of Grassmann 
variables and a minimal basis of composite fields has been set up which reproduces 
the full pairing spectrum in the degenerate case together with the associated wave 
functions. 

When many single particle levels are present, new facets of the pairing physics 
disclose themselves. 

First a third type of excitations (promotion of pairs to higher levels) can occur. 
These have been classified in terms of a new quantum number, the "like-seniority" , 
closely related to the Gaudin number j^HI) which has been shown to coincide with the 
true seniority when the strength of the interaction is large. Second, critical values 
of the coupling constant G appear, signalling the transition between two different 
regimes: one dominated by the mean field, the other by the pairing bosonizing inter- 
action. Notably the number of the critical values of G, which is crucially linked to 
the exsistence of a lowest single particle level with pair degeneracy Q = 1, increases 
with the number of pairs: however all the Gcr collapse into a common value when 
the number of levels tends to infinity. 

The next step in the study of the pairing Hamiltonian will be to consider the 
general case of n pairs living in a set of single particle levels and explore whether 
the signatures of a Goldstone boson and of the Higgs excitation are still present, so 
shading further light on the phase transition to superconductivity occurring in finite 
and infinite Fermi systems. In this connection the escape mechanism of the pairs 
from the grid of the single particle levels should be carefully examined in relation 
to the pair degeneracy of the latter. 

Finally, worth to be explored are new scheme recently proposed for the bosoniza- 
tion of finite Fermi systems, specifically of atomic nuclei: in particular the one 
outlined in Refs. j;-i4| which, always in the framework of the path integral, in- 
troduces the concept of boson dominance in computing the trace expressing the 
partition function of a nucleus. 

Acknowledgements 

Most of the work presented here is the fruit of a longstanding collaboration with 
R. Cenni, A. Molinari and F. Palumbo, whom we would like to thank for many 
useful discussions and careful reading of the manuscript. 

References 

[1] J. Bardeen, L. N. Cooper and J. R. Schrieffer, Phys. Rev. 106, 162 (1957). 

[2] A. Bohr, B. R. Mottelson, D. Pines, Phys. Rev. 110, 936 (1958). 

[3] A. Arima, T. Otsuka, F. lachello and I. Talmi, Phys. Lett. , B 66, 201, 1977; 
T. Otsuka, A. Arima, F. lachello and I. Talmi, Phys. Lett, B 76:139, 1978. 



42 M.B. Barbara, M.R. Quaglia 



[4] J. Bardeen, L. N. Cooper and J. R. SchriefFer, Phys. Rev. 108 (1957) 1175. 

[5] J.N. Ginocchio Phys. Lett. B 79, 173, 1978; J.N. Ginocchio Nucl. Phys. A 
376, 438, 1982; H.B. Geyer Phys. Rev. C 34, 2373, 1986; C.W. Johnson and 
J.N. Ginocchio Phys. Rev. C 50, 571, 1994. 

[6] Algebraic Approaches to Nuclear Structure, edited by R.F. Casten, Harwood 
Academic Pubhshers, 1993. 

[7] S. Weinberg, The quantum theory of fields, Vol.2, Cambridge University Press, 
1996. 

[8] M. C. Cambiaggio, A. M. F. Rivas and M. Saraceno, larXiv:nucl-th/970803Tl 

[9] A. K. Kerman, R. D. Lawson and M. H. Macfarlane, Phys. Rev. 124, 162 
(1961). 

[10] A. K. Kerman, Ann. Phys. 12, 300 (1961). 
[11] P. W. Anderson, Phys. Rev. 112, 1900 (1958). 

[12] M. B. Barbaro, A. Mohnari, F. Palumbo and M. R. QuagHa, Phys. Rev. C 70, 
034309 (2004). 

[13] F. Palumbo, Phys. Lett. B328 (1994) 79 ; G. De Franceschi and F. Palumbo, 
Modern Phys. Lett. All (1995) 901. 

[14] M. B. Barbaro, A. Mohnari and F. Palumbo, Nucl. Phys. B 487, 492 (1997) 
|arXiv:hep-th/960716Ql . 

[15] P. W. Anderson, Phys. Rev. 130, 439 (1963). 

[16] S. Fubini and A. Molinari, Annals Phys. 221 (1993) 1. 

[17] M. B. Barbaro, A. Molinari, F. Palumbo and M. R. Quagha, Phys. Lett. B 
476, 477 (2000) arXiv:nucl-th/9912035 . 

[18] M. B. Barbaro, A. Molinari, F. Palumbo and M. R. Quagha, Atti Accad. Sci. 
Torino. Sci. Fis. Mat. Natur. 133, 3 (1999) arXiv:nucl-th/0005055 . 

[19] P. Ring and P. Shuck, The Nuclear Many Body Problem, (Springer- Verlag, 
New York, 1980), p. 225. 

[20] A. Bohr and B.R. Mottelson, Nuclear Structure, Vol. II, (Benjamin, New York, 
1975). 

[21] G. Sierra, J. Dukelsky, G. G. Dussel, J. von Delft and F. Braun, Phys. Rev. B 
61, 11890 (2000) arXiv:cond-mat /99090 1 5 . . 



Bosonization of the Pairing Hamiltonian 



43 



[22] J.J3ukelsky, C. Esebbag and P. Schuck, Phys. Rev. Lett. 87 (2001) 066403 
|aIXiv :cond-mat/01 07477| . 

[23] R.W. Richardson and N. Sherman, Nucl. Phys. , 52, 221, 1964; R.W. Richard- 
son, Phys. Rev. , 141,949, 1966. 

[24] M. B. Barbaro, L. Fortunate, A. Molinari and M. R. QuagHa, Phys. Rev. 
C 64, 011302 (2001) arXiv:nucl-th /0105006 ; M. B. Barbaro, L. Fortunato, 
A. Mohnari and M. R. Quagha, arXiv:nucl-th/0012004, 

[25] M. B. Barbaro, R. Cenni, A. Molinari and M. R. Quagha, Eur. Phys. J. A 22, 
377 (2004) ,arXiv:nuc l-th/0306062,. 

[26] M. Gaudin, "Modeles Exactement Resolus", les Editions de Physique, France, 
1995. 

[27] J. M. Roman, G. Sierra and J. Dukelsky, Phys. Rev. B 67 (2003) 064510. 

[28] M. Barbaro, R. Cenni, A. Mohnari and M. R. Quagha, Phys. Rev. C 66, 034310 
(2002) arXiv:nucl-th/0206070j. 

[29] D.J. Rowe, Nuclear Collective Motion, Methuen and Co. LTD, London, 1970. 

[30] G. Sansone and J. Gerretsen, Lectures on the theory of functions of a complex 
variable, P. Noordhoff, 1960. 

[31] J. M. Roman, G. Sierra and J. Dukelsky, Nucl. Phys. B 634, 483 (2002) 
|arXiv:cond-iSS;/ 0202070| . 

[32] E. A. Yuzbashyan, A. A. Baytin and L. Altshuler, Phys. Rev. B 68 (2003) 
214509. 

[33] M. Hasegawa and S. Takaki, Phys. Rev. C 35 (1987) 1508. 
[34] F. Palumbo, arXiv:nucl-th/0412042| 
[35] F. Palumbo, |arXiv:nucl-th /0405045| 



